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A procedure is developed for using nonlinear experimental response data to guide the 
modal basis selection in a nonlinear reduced-order simulation. The procedure entails using 
nonlinear acceleration response data to first identify proper orthogonal modes. Special 
consideration is given to cases in which some of the desired response data is unavailable. 
Bases consisting of linear normal modes are then selected to best represent the 
experimentally determined transverse proper orthogonal modes and either experimentally 
determined in-plane proper orthogonal modes or the special case of numerically computed 
in-plane companions. The bases are subsequently used in nonlinear modal reduction and 
dynamic response simulations. The experimental data used in this work is simulated to 
allow some practical considerations, such as the availability of in-plane response data and 
non-idealized test conditions, to be explored. Comparisons of the nonlinear reduced-order 
simulations are made with the surrogate experimental data to demonstrate the effectiveness 
of the approach. 


C, C 

f, f 

3 > 

/. 

M, M 

v 

CO 

p 

R x 

x, q 
X 

x 


Nomenclature 

system damping matrix in physical and modal degrees of freedom 
force excitation vector in physical and modal degrees of freedom 
nonlinear restoring force vector in physical and modal degrees of freedom 

modal basis matrix 

proper orthogonal value (POV) matrix 

system mass matrix in physical and modal degrees of freedom 
cumulative modal amplitude participation factor 
undamped natural frequency 
proper orthogonal mode (POM) matrix 
displacement correlation matrix 

displacement response vector in physical and modal degrees of freedom 
displacement snapshot matrix 
modal amplitude participation factor 
viscous damping factor 


I. Introduction 

N ONLINEAR random response analysis of large structural components in physical degrees of freedom (DoFs) 
can be associated with a prohibitive computational cost. Reduced-order modeling has been shown to be a 
viable alternative for a range of loading and response conditions. 1 " 8 For a typical aerospace structure having 
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complex geometry and loading conditions (combined mechanical, thermal, acoustic, and pressure), selection of basis 
functions through which the reduced-order system is formed is usually not intuitive. A rigorous modal basis 
selection procedure resulting in a computationally efficient and accurate reduced-order simulation was recently 
developed/’ As part of that procedure, system identification was performed and a criterion for modal basis selection 
was established. In that work, the system identification procedure used nonlinear dynamic response data from a full- 
order numerical simulation in physical DoFs. The data consisted of both transverse and in-plane DoFs. The goal of 
this work is to explore the use of experimental data in the system identification procedure. A key benefit of the 
approach will be the elimination of a full-order simulation currently required for the system identification procedure. 

An aluminum beam structure under a transverse, point-load, random excitation is used to demonstrate the 
approach. Two variations are considered: an idealized case with perfectly clamped boundary conditions and known 
transducer mass and location, and a non-idealized case with unknown compliant boundary conditions and perturbed 
transducer mass and locations. For the idealized case, nonlinear acceleration response data is used in a proper 
orthogonal decomposition (POD) analysis to obtain proper orthogonal modes (POMs). Since in-plane acceleration 
data is not typically measured, transverse acceleration data are used to identify only the transverse POMs. These, in 
turn, are used as static displacement fields in a nonlinear static analysis to determine corresponding in-plane 
behavior. A load-independent modal basis is subsequently identified consisting of a set of linear normal modes 
(LNMs) resembling the combination of experimentally determined transverse POMs and numerically determined in- 
plane displacement shapes. An alternative approach, using experimentally determined in-plane POMs instead of 
numerically determined displacement shapes, is also demonstrated. For the non-idealized case, the linear dynamics 
model is first updated to match the linear vibration response. The model update procedure uses linear static 
deflection data to adjust the stiffness, and the fundamental frequency to adjust the mass. The basis selection 
considered uses experimentally determined transverse POMs and numerically determined in-plane displacement 
shapes. Note that the nonlinear modal reduction used by the authors does not statically condense in-plane DoFs. 3 
Therefore, a technique similar to the current approach, 10 but which identifies only transverse LNMs, would generate 
only part of the required basis. 

Because the focus of this work is on methods development, simulated experimental data is used in lieu of actual 
experimental data. This surrogate experimental data offers the ability to explore practical considerations including 
the availability of in-plane acceleration data and the effect of modeling uncertainty on the simulated response. 
Further, it serves to guide further experimental studies insofar as measurement specifications are concerned. 

II. Nonlinear Modal Reduction and Simulation 

Since the modal reduction utilizing the indirect approach employing a nonlinear stiffness evaluation procedure 
was detailed by the authors in several recent publications, 1 ' 1 1-14 only an abbreviated description is presented herein. 
The equations of motion of the nonlinear full-order system in physical DoFs may be written as 

Mx(t) + Cx0) + f N| (x(0) = f(0 (1) 

where M and C are the system mass and damping matrices, x is the displacement response vector and f is the force 
excitation vector, respectively. The nonlinear restoring force f NL is a vector function, which generally includes the 
linear, quadratic, and cubic stiffness terms. By applying the modal coordinate transformation 

x(0 = Oq(t) (2) 

to Eq. (1), a modal equation of motion can be written as 

Mq(0 + Cq(t) + f NL (q,(t), q , (0, • • • , q, (0) = f« (3) 

where q=[^ 1 ,^ 2 ,...,^ i ] T is a generalized coordinate vector and ® is the set of L selected LNMs. For mass- 
normalized eigenvectors, 

M = 0> T M® = riJ c = <5 t CO = J (4) 
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and co are the undamped natural frequencies and C r are the viscous damping factors. The modal excitation force is 

f = <t'f . 

Since the nonlinear restoring force f is generally not known in the context of a commercial finite element (FE) 

program, an indirect means of evaluating the nonlinear stiffness is required. Different combinations of scaled LNMs 
can be used to form a set of prescribed displacement fields. Using a nonlinear FE static analysis, nonlinear restoring 

forces f can be computed in physical DoFs and transformed to the generalized coordinates as f NL = O I M . The r th 

nonlinear modal restoring force can be written as a summation of modal linear, quadratic, and cubic stiffness 
components as 


L L L L L L 

f r NL (?, , q 2 , ■ ■ ■ ,q L ) = Z + Z Z + Z Z Z r = \, 2,..., l. (5) 

j = 1 j=l k=j j = 1 k=j l=k 

Since the left-hand-side of the Eq. (5) is known, as are the assumed scaling factors q on its right-hand-side, solution 
of a simple algebraic system of equations is only needed to arrive at the linear d , quadratic a, and cubic b modal 
stiffness coefficients. Eq. (3) is numerically integrated using a 4 ,h -order Runge-Kutta scheme and physical 
displacements are obtained through the inverse modal coordinate transform. 

III. System Identification and Modal Basis Selection Procedure 

Three approaches were previously offered for selecting a modal basis following the system identification 
procedure. 9 These were the modal amplitude participation (MAP), modal energy participation (MEP), and estimated 
POM frequency (EPF) procedures. Since the present study will consider a hybrid of experimentally determined 
transverse POMs and numerically determined corresponding in-plane behavior, the MAP procedure, which 
considers these DoF types separately, will be utilized. 

When physical DoFs are chosen to characterize a response, a data set can be formed as an accumulation of n 
instantaneous displacement fields to produce a displacement snapshot matrix X . The snapshot matrix contains the 
selected set of N DoFs resulting in its size being n x N. The sample rate and spatial resolution of the snapshot 
matrix must be sufficient to resolve the system’s temporal and spatial characteristics of interest. The displacement 
correlation matrix R x of size N x N is formed as 

R x =-x‘x (6) 

n 


An eigenanalysis of the displacement correlation matrix R x is next performed, i.e., 

[R x -M]p = 0 (7) 

to obtain the POM matrix P = [p p, ... p ] and the diagonal proper orthogonal value (POV) matrix, /. , both of 

size N x N. 15 Note that an alternative means of obtaining similar quantities is via singular value decomposition. 1 6 
The contribution of each POM to the overall dynamic response is given by 

/ i = l,...,N (8) 

where % is the f' POM modal amplitude participation factor. The sum of all POM modal amplitude participation 
factors is unity. When the dominant M POMs are selected, their cumulative participation, v , can be expressed as 
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0 < u < 1 . 


(9) 


M 

jM 

Retention of only M selected POMs reduces the size of P to N x M. The MAP approach requires as many criteria as 
there are DoF types used in the system identification, e.g. one cutoff value for the transverse displacement and one 
value for the in-plane displacement. Using such cutoff values, a set of POMs may be selected based on the modal 
amplitude participation factor, the cumulative participation, or some combination thereof. 

The direct use of POMs as the basis functions is not preferred as they are specific to the loading under which 
they were determined and, consequently, applicability of such a basis to alternative loading conditions may be 
limited. Instead, a more robust basis, formed from load-independent LNMs <l> which resemble the M selected 
POMs, is sought. The authors previously identified two means of doing so; 11 one using the modal assurance 
criterion (MAC) and the other using a modal expansion. ' ’ The MAC approach was used in this work. Fora 
pair of POM p and normal mode (p vectors, the MAC value can be obtained from 


MAC(p t , <p ( ) 



T 

P 

<P, 

2 

(pj 

Pj 

9- 

t- 1 

9- 


k - 1, 

/ = 1 , 


M 
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Once the MAC matrices are determined, a MAC threshold level is used to identify the normal modes best 
resembling the selected POMs. A threshold of MAC > 0.5 was used in this work. The end result of the basis 
selection process is a set of L selected LNMs, L«N, for use in the nonlinear modal reduction. 

IV. Idealized Case 

The idealized case was used to develop the procedure for generating transverse and in-plane POMs from 
acceleration data, and for generating numerically determined in-plane POM replacements in the absence of in-plane 
acceleration data. A planar aluminum beam structure, measuring 0.4572 m in length with a cross-section of 
25.4 mm wide by 2.286 mm high, was previously investigated by the authors 13 and was the subject of the current 
study. The beam was perfectly clamped at both ends and subject to transverse point-load random excitation at the 
14- span location, as shown in Figure 1. The zero-mean Gaussian distributed excitation had a flat spectrum to 
1500 Hz. Two loading levels were considered; a low level loading of 0.1N root-mean-square (RMS) was used to 
excite the beam in the linear response regime, and a high level loading of 12. 8N RMS was used to generate a 
geometrically nonlinear response. 


38.1 mm 



A. Surrogate Experiment 

As previously indicated, the “experiments” performed in this work were actually full-order (physical DoF) 
numerical simulations of the dynamic response under virtual test conditions. The full-order finite element model 
was comprised of 144 ABAQUS ls two node B21 elements, each measuring 3.175 mm in length. The B21 element 
features three DoFs per node; one rotational DoF and two translational (transverse and in-plane) DoFs. Material 
properties were: Young’s modulus E = 73.1 1 GPa, shear modulus G = 27.59 GPa, and mass density p = 2763 kg/m 3 . 
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Mass proportional damping was specified as 14.54 s' 1 , which corresponded to a viscous damping factor of 2.78% at 
the linear fundamental frequency of 41.65 Hz. 

Accelerometers were modeled as 0.4g lumped masses, and were uniformly spaced 38.1 mm apart, as shown in 
Figure 1. Additionally, a 25g load cell was modeled as a lumped mass at the 14-span drive point. The stinger and 
shaker shown in the schematic were not modeled. Acceleration time histories were simultaneously recorded at a 
sampling rate of 20k samples/s. Experimental data was simulated for 104 s, with the first 0.5 s discarded to remove 
the startup transient. Transverse and in-plane acceleration power spectral densities (PSDs) at the 14-span drive point 
are shown in Figure 2 and Figure 3, respectively. A 32,768-point transform was used in these and subsequent PSDs, 
giving a frequency resolution of 0.61 Hz. The nonlinear response is evident in the broad peaks and significant 
response outside the 0-1500 Hz excitation bandwidth. Note that the in-plane acceleration response rolls off at twice 
the excitation bandwidth (not shown) due to frequency doubling. Frequency doubling occurs because each half 
cycle of upward and downward transverse motion results in a frill cycle of in-plane stretching. 




Figure 2: Nonlinear transverse drive point Figure 3: Nonlinear in-plane drive point 

acceleration PSD for the 12. 8N RMS loading. acceleration PSD for the 12. 8N RMS loading. 


While the number of active DoFs in the finite element model is 429 (143 nodes x 3 DoF/node), the number of 
active DoFs in the experiment was either 11 when using uni-axial accelerometers to “measure” transverse 
acceleration, or 22 when using bi-axial accelerometers to “measure” both transverse and in-plane acceleration. A 
practical consideration in the design of the experiment therefore is that the number of measurements must be greater 
than or equal to the anticipated number of M selected POMs. 


B. System Identification and Modal Basis Selection 

The POD analysis, as discussed in Section III, uses displacement time histories to form the correlation matrix, 
R x , not the acceleration time histories. Therefore, the first step in the analysis was to integrate the nonlinear 

acceleration time history, obtained from the 12.8N RMS loading, to obtain the displacement time history. Shorter 
records of duration 1.6384 s were used to form the 32,768-point («) displacement snapshots used in the POD 
analysis. The integration operation was performed in the frequency domain by taking the complex 32,768-point 

acceleration transform, dividing by -co 1 , mirroring the complex conjugate about the Nyquist frequency, and taking 
the inverse transform to obtain the displacement. Leakage was avoided by selecting an analysis window equal to the 
record length, hence windowing was not required. The resulting displacement time history was then high-pass 
filtered with a 3 rd -order Butterworth filter to remove the DC component. The MATLAB 19 filter function “filtfilt” 
was used to eliminate phase distortion. The resulting transverse and in-plane displacement response PSDs are 
shown in Figure 4 and Figure 5, respectively. 
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Figure 4: Nonlinear transverse drive point Figure 5: Nonlinear in-plane drive point 

displacement PSD for the 12. 8N RMS loading. displacement PSD for the 12. 8N RMS loading. 

1. Transverse Modes 

In this study, transverse POMs were always determined from integrated transverse accelerometer data. These 
POMs served two purposes. The primary purpose was to select part of the LNM basis using the MAC analysis 
discussed in Section III. The second purpose was to serve as a set of prescribed transverse deflections in a nonlinear 
static analysis to numerically determine the corresponding in-plane displacement fields. These displacement fields 
act as alternatives to in-plane POMs when in-plane accelerometer data is unavailable. Because of the requirement to 
prescribe the transverse deflections at all the nodes and not just at the accelerometer locations, the transverse POMs 
were first fit with a cubic spline. Thus, another practical consideration is that an appropriate spatial distribution of 
accelerometer locations is required to characterize the POM shape. The spline-fit would obviously not work well, 
for example, if all the accelerometer locations were toward one side. For consistency, the fitting operation was 
performed regardless of whether or not in-plane accelerometer data was available. As such, the subsequent MAC 
analysis was performed using all 145 transverse DoFs ( N ). Had the requirement to derive an in-plane POM 
alternative not been present, the LNM shapes obtained from the eigenanalysis of the full-order model could have 
been decimated to the 1 1 accelerometer locations. 

The first 6 POMs, ranked according to their participation, as computed per Eq. (8), are shown in Figure 6 - 
Figure 11. In each figure, the raw experimentally determined POM is shown at the nodal spacing corresponding to 
the accelerometer locations. The spline-fit POM used for both the MAC and nonlinear static analyses are shown at 
each nodal location. For comparison, POMs obtained directly from the full-order FE analysis displacement time 
histories are also shown. These POMs are included to contrast the present experimental approach with the former 
simulation-based approach. 9 ' 1 1-13 It is seen that the set of POMs from the two methods compare well. 




Figure 7: Transverse POM 2. 
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Node Number 

Figure 8: Transverse POM 3. 



Node Number 

Figure 10: Transverse POM 5. 



Figure 9: Transverse POM 4. 



Figure 11: Transverse POM 6. 


System identification data for the 11 transverse POMs obtained from acceleration data is presented in Table 1. 
Shown are the POV participation factors computed per Eq. (8), the cumulative POV participation factors computed 
per Eq. (9), the corresponding LNM (computed from the eigenanalysis of the full-order model) as determined by the 
MAC analysis, the LNM natural frequency, and the MAC values computed per Eq. (10). The LNM number was 
assigned in order of ascending natural frequency. It is seen that inaccuracies in the experimentally obtained POMs 
did not result in a low MAC value, nor by a need to represent a single POM with multiple LNMs. 


Table 1: System identification and basis selection using transverse POVs obtained from acceleration. 


POM 

POV (%) 

Cum POV (%) 

LNM 

LNM Natural Freq (Hz) 

MAC 

1 

79.20 

79.20 

1 

41.647 

0.997 

2 

16.37 

95.58 

2 

104.78 

0.919 

3 

3.35 

98.92 

3 

219.14 

0.862 

4 

0.61 

99.53 

6 

762.42 

0.763 

5 

0.32 

99.85 

4 

394.57 

0.986 

6 

0.086 

99.94 

5 

587.94 

0.941 

7 

0.056 

99.99 

7 

1041.8 

0.590 

8 

3.7e-4 

99.99 

8 

1395.9 

0.910 

9 

7.4e-5 

99.99 

9 

1701.8 

0.725 

10 

1.9e-5 

99.99 

10 

1919.2 

0.770 

11 

4.0e-8 

99.99 

11 

2244.1 

0.865 


The selection of the transverse portion of the modal basis for the reduced order analysis was therefore based 
exclusively on identifying POMs above some threshold value of the cumulative POV factor. A value of v >99.99% 
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was chosen for this purpose resulting in selection of 7 POMs (M r ) and 7 corresponding LNMs (L r ) for the transverse 
portion of the modal basis, indicated by the red-line in Table 1. The 7 selected POMs are also used to numerically 
determine corresponding in-plane behavior in the absence of in-plane acceleration data. Finally, it should be noted 
that the POD/MAC analysis using finite element displacement data identified the same 7 transverse LNMs. 

2. In-Plane Modes 

As previously indicated, two methods for selecting the in-plane LNM basis are employed. The first relies on in- 
plane acceleration, as would be acquired with transverse acceleration using a biaxial accelerometer. In this case, the 
processing of in-plane acceleration to displacement time histories, and the POM/MAC analyses are identical to the 
approach presented in the previous section. In the more likely scenario, in-plane acceleration data is not available 
and the in-plane response must be derived from the transverse POMs. Both methods are next discussed. 

For the case when in-plane acceleration data is available, the first 6 POMs are shown in Figure 12 - Figure 17. 
In each figure, the raw experimentally determined POM is shown at the nodal spacing corresponding to the 
accelerometer locations. The spline-fit POMs used for the MAC analyses are shown at each nodal location. For 
comparison, POMs obtained directly from the full-order FE analysis displacement time histories are also shown. It 
is seen that the set of POMs from the two methods are similar, although agreement is not as good as was observed in 
the transverse POMs (see Figure 6 - Figure 1 1). 



Node Number 


Figure 12: In-plane POM 1. 





Figure 14: In-plane POM 3. 


Figure 15: In-plane POM 4. 
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Node Number Node Number 

Figure 16: In-plane POM 5. Figure 17: In-plane POM 6. 

System identification data for the in-plane POMs obtained from acceleration data is presented in Table 2. It is 
seen that inaccuracies in the experimentally obtained POMs did result in lower MAC values compared to the 
transverse POMs. In nearly half the cases, the MAC criterion was unable to identify a LNM using the same cutoff 
value of 0.5. As such, the cumulative POV used only reached 97.64% for the 6 LNMs (Z/) identified. For this case, 
use of an alternative mapping technique, e.g. modal expansion, 9 * ' 11 ’ 12 might have identified a greater number of 
LNMs and hence higher cumulative POV than the MAC criterion employed. While the MAC analysis using POMs 
obtained directly from the finite element analysis displacement time histories (not shown) was able to identify three 
additional LNMs (34, 39 and 69), it could not identify LNMs corresponding to POMs 9 and 10. 


Table 2: System identification and basis selection using in-plane POVs obtained from acceleration. 


POM 

POV (%) 

Cum POV (%) 

LNM 

LNM Natural Freq (Hz) 

MAC 

1 

79.00 

79.00 

13 

3914.9 

0.953 

2 

9.11 

88.10 

20 

7411.1 

0.796 

3 

5.61 

93.71 

31 

17354 

0.508 

4 

3.41 

97.12 

26 

12289 

0.642 

5 

1.57 

— 

— 

— 

<0.5 

6 

0.65 

— 

— 

— 

<0.5 

7 

0.29 

97.41 

45 

33492 

0.658 

8 

0.23 

97.64 

42 

28671 

0.654 

9 

0.08 

— 

— 

— 

<0.5 

10 

0.05 

— 

— 

— 

<0.5 

11 

0.01 

— 

— 

— 

<0.5 


For the more typical case when in-plane acceleration data is unavailable, the selected set of 7 experimentally 
determined transverse POMs were individually applied as static displacement fields in a series of nonlinear static FE 
analyses. The POMs were scaled to a maximum displacement of 5 mm, or about two times the beam thickness, as 
this was close to the maximum transverse displacement found from the integrated acceleration data. Six of seven of 
the resulting displacement fields, heretofore referred to as companions, are shown in Figure 18 - Figure 23. It is 
clear that the companion shapes are vastly different than the POM shapes found from in-plane acceleration data. 
However, a linear superposition of companion modes can yield shapes similar to the in-plane POMs. It is seen that 
superposition of companions 1-7 reconstruct shapes similar to in-plane POMs 1 and 2, as shown in Figure 24 and 
Figure 25, respectively. Recall that these two POMs account for 88% of the in-plane response. Therefore, it is 
expected that similarities will exist between the LNMs identified using companions and those identified using the 
in-plane POMs. 

Basis selection data for the in-plane companions is presented in Table 3. Like the in-plane POMs, MAC 
analyses using the in-plane companions resulted in lower MAC values compared to the transverse POMs. Because 
the companions were not obtained from a POD analysis, the individual and cumulative POV participation factors are 
not available. Therefore, the cutoff was effectively established by the cutoff used to select the 7 transverse POMs. 


9 

American Institute of Aeronautics and Astronautics 


A consequence of using companions instead of POMs in the MAC analyses was the identification of only three 
unique in-plane LNMs. However, these three LNMs are a subset of those identified using in-plane POMs, and are a 
reflection of the fact that companions can be superposed to reconstruct shapes similar to in-plane POMs. The effect 
of using only three in-plane LNMs (L f ) on the reduced-order simulation results, versus 6 and 9 for the acceleration 
and displacement-based POMs, respectively, will next be discussed. 




Figure 18: In-plane companion 1. 



Figure 20: In-plane companion 3. 



Figure 19: In-plane companion 2. 




Figure 22: In-plane companion 5. 


Figure 23: In-plane companion 6. 
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Figure 24: Superposition of companions for Figure 25: Superposition of companions for 

reconstruction of in-plane POM 1. reconstruction of in-plane POM 2. 


Table 3: Basis selection using in-plane companions obtained from transverse POMs. 


Companion 

LNM 

LNM Natural Freq (Hz) 

MAC 

1 

31 

17354 

0.562 

2 

20 

7411.1 

0.619 

3 

20 

7411.1 

0.557 

4 

13 

3914.9 

0.946 

5 

— 

— 

<0.5 

6 

20 

7411.1 

0.842 

7 

13 

3914.9 

0.901 


C. Reduced-Order Simulation Results 

Four different bases were considered for the reduced-order analyses. The first consisted of only 7 transverse 
LNMs derived from the transverse acceleration data, and is subsequently referred to as the 7-mode acceleration 
derived basis. As will be shown, this basis was included to show the importance of in-plane modes and not as a 
desired basis. The second basis consisted of 7 transverse LNMs and 6 in-plane LNMs, each derived from 
acceleration data. This is subsequently referred to as the 13 -mode acceleration derived basis. The third consisted of 
7 transverse LNMs derived from acceleration data, and 3 in-plane LNMs derived from companions. This is 
subsequently referred to as the 10-mode hybrid basis. The last basis consisted of 7 transverse LNMs and 9 in-plane 
LNMs, each derived from displacement data. This is subsequently referred to as the 16-mode displacement derived 
basis. Its inclusion allows comparisons of the new acceleration derived and hybrid bases, with the previous 
approach of using numerical simulation data for the basis selection. Table 4 gives a breakdown of selected 
transverse and in-plane LNMs for each basis. 

Table 4: Content of each modal basis. 


Basis 

Transverse LNMs 

In-Plane LNMs 

7 mode acceleration derived 

1-7 

none 

10 mode hybrid basis 

1-7 

13,20,31 

13 mode acceleration derived 

1-7 

13,20,26,31,42,45 

16 mode displacement derived 

1-7 

13,20,26,31,34, 39, 42,45, 69 


A comparison of the transverse displacement PSDs at the drive point is shown in Figure 26. All reduced-order 
analyses compare well with the simulated experimental data, with differences noted at the higher frequencies due to 
modal truncation. Because of nonlinear coupling, peaks in the response are seen outside of the excitation 
bandwidth. To highlight the differences between the reduced-order analyses, a close-up of the low frequency 
response is shown in Figure 27. There it is seen that a lesser number of in-plane modes in the basis (recall that all 
reduced-order analyses have the same 7 transverse modes) acts to over-predict the degree of nonlinear spring- 
hardening. The 7-mode acceleration derived basis, with no in-plane modes, shows the most spring-hardening as 
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indicated by a shifting of peaks to higher frequencies, a greater degree of peak broadening, and reduced amplitude 
relative to the experiment. Addition of in-plane modes lessens the degree of nonlinear spring-hardening, with the 
13 -mode acceleration derived basis and 16-mode displacement derived basis virtually indistinguishable from the 
experiment. 

The in-plane displacement PSDs at the drive point are shown in Figure 28, and in close-up view in Figure 29. 
All three reduced-order analyses compare favorably with the experiment at low frequencies and follow similar 
trends at higher frequencies. It should be noted that the behavior shown is solely due to coupling between the 
transverse and in-plane modes, not direct excitation of the in-plane modes which are at higher frequencies (see Table 
2 and Table 3). The 7-mode acceleration derived basis is incapable of predicting any in-plane response since it does 
not contain any in-plane modes. The close-up view indicates that the reduced-order analyses improve by increasing 
the number of modes, as expected. 



Figure 26: Comparison of transverse drive point 
displacement PSD for the 12. 8N RMS loading. 



Figure 28: Comparison of in-plane drive point 
displacement PSD for the 12. 8N RMS loading. 



Figure 27: Close-up of transverse drive point 
displacement PSD for the 12. 8N RMS loading. 



Frequency (Hz) 


Figure 29: Close-up of in-plane drive point 
displacement PSD for the 12. 8N RMS loading. 


Finally, recall that the basis consisting of load-independent LNMs, instead of one consisting of load-dependent 
POMs, was selected in order to form a robust basis applicable over a range of loadings. This is demonstrated in 
Figure 30 and Figure 31 for a far lower excitation of 0.1N RMS, where excellent comparison of all reduced-order 
analyses with the simulated experimental data is seen across the frequency range. Again, exceptions are noted at 
high frequency due to modal truncation, and in the 7-mode basis for the in-plane response. 
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The conclusion for this idealized case is that reduced-order models formed with LNM bases identified from a 
system identification approach using either both transverse and in-plane acceleration data, or only transverse 
acceleration data, are a viable alternative to bases identified with the previous displacement based approach. Of the 
two approaches, the basis selected using both transverse and in-plane acceleration data was more accurate as it was 
able to identify a greater number of in-plane modes. Next investigated is the non-idealized case where modeling 
uncertainty in the boundary conditions and accelerometer mass loading are considered to more closely mimic actual 
experimental conditions. 




Figure 30: Comparison of transverse drive point Figure 3 1 : Comparison of in-plane drive point 

displacement PSD for the 0. IN RMS loading. displacement PSD for the 0. IN RMS loading. 

V. Non-idealized Case 

The non-idealized case utilized the same beam structure as the idealized case, but with compliant boundary 
conditions and variations in the accelerometer mass distribution. The same transverse random loading was applied 
at the 14-span location. System identification was performed at the 12. 8N RMS loading level. The basis used in the 
reduced-order model was developed from the system identification using transverse acceleration data alone, i.e. the 
hybrid basis, since in-plane acceleration data is not typically available. 

A. Surrogate Experiment 

Compliant boundary conditions were simulated using two node CONN2D2 connector elements in ABAQUS. 18 
The in-plane and transverse spring stiffness were set high at lOGN/m to simulated immovable boundaries. The 
rotational stiffness was set at 500 N-m on the left hand side, and 100 N-m on the right hand side. 

The mass loading of the accelerometers was altered in both magnitude and distribution. The magnitude was 
changed to account for varying lead wire mass by allowing the total mass of each accelerometer to vary between 
0.4-0.6g. This increased the total added mass from 44g in the idealized case to 55.62g for the non-idealized case. 
The accelerometers were assumed to have a mounted base dimension of 5.8 mm. They therefore span two elements 
or 3 nodes in the surrogate experiment finite element model. The randomized mass was itself randomly distributed 
across the 3 nodes to allow the center mounting location to deviate from the nominal location by as much as 1.6 mm. 
The load cell was assumed to have a mounted base dimension of 16 mm and its mass was distributed equally across 
5 nodes. The total mass of 25g was left unchanged from the idealized case. 

The net effect of these changes was a 12% reduction in the fundamental frequency from 41.65 Hz to 36.72 Hz, 
and asymmetry of the LNM shapes (not shown) due to different compliance at the left and right boundaries. Note 
that the compliant boundary conditions result in a fundamental frequency much closer to an idealized clamped 
condition than to an idealized simply supported condition (17.22 Hz). 

B. Full-Order Model Calibration 

Before embarking on the reduced-order model, it was first necessary to calibrate the full-order model from which 
the ROM is derived. This was not necessary for the idealized case because the full-order model was the same as the 
surrogate experimental FE model. Here, the intention is to calibrate the full-order model in a manner that might be 
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performed by the experimental analyst. In this case, the calibration method is focused on matching the fundamental 
frequency, without regard to the mode shape. In this process, material properties of the beam were assumed to be 
well known and were not altered. 

A static deflection test was first performed to calibrate the boundary stiffness. Via simulated experiment, a static 
transverse load of 1.0 N was applied to the surrogate experimental FE model at the !4-span shaker location and a 
transverse deflection of 0.126 mm was read using a linear variable differential transformer (LVDT). The rotational 
stiffness of the CONN2D2 elements in the full-order model were adjusted to equal values until the simulated LVDT 
measurement matched that of the experiment. No attempt was made to use mode shapes to make a differential 
adjustment. Ultimately, a rotational stiffness of 361.94 N-m was equally applied to the left and right sides. 

A dynamic test was next performed to calibrate the effect of mass. In the full-order model, all masses were 
applied in an idealized fashion, i.e. at a single node. Only the mass of the accelerometers was adjusted to match the 
fundamental frequency, and with the same adjustment applied to all masses. It was determined that an 
accelerometer mass of 0.61g was required to match the fundamental frequency. 

The accelerance frequency response functions of the experiment and calibrated full-order models are shown at 
the drive point and center span locations in Figure 32 and Figure 33, respectively. The calibrated full-order model 
matches the experimental data very well below 300 Hz. At higher frequencies, the softness of the full-order model 
relative to the experiment becomes more apparent. Since the reduced-order model is derived from the calibrated 
full-order model, similar behavior is expected. 




Figure 32: Comparison of transverse accelerance at Figure 33: Comparison of transverse accelerance at 

drive point. center s P an - 

C. System Identification and Modal Basis Selection 

The hybrid approach using transverse acceleration data and companions for the in-plane response yielded a 
modal basis consisting of 7 transverse modes and 2 in-plane modes, as indicated in Table 5. Like the results in 
Table 3 for the idealized case, the same LNMs were identified with multiple companions and the MAC values were 
generally low. 

D. Reduced-Order Simulation Results 

A comparison of results from the calibrated reduced-order model and the experimental are shown in Figure 34 - 
Figure 37. For reference, also shown are the simulated responses from the calibrated full-order model. These 
represent the best any reduced-order model can be expected to perform. As would be expected from the linear 
accelerance in Figure 32, the full-order and reduced-order analyses compare most favorably at low frequencies and 
differ with increasing frequency. From the close-ups in Figure 35 and Figure 37, the calibrated full-order model is 
seen to compare very well with the experiment. Like the idealized case, the reduced order model indicates 
somewhat greater nonlinearity, with a quality comparable to that of the 10-mode basis for the idealized case. 
Therefore, it may be concluded that modeling uncertainties do not affect the reduced-order model results any more 
than the full-order model from which it was derived. 
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Table 5: System identification and basis selection results using the hybrid approach. 


POM/Companion 

POV (%) 

Cum POV (%) 

LNM 

LNM Natural Freq (Hz) 

MAC 

POM 1 

72.18 

72.18 

1 

36.721 

0.980 

POM 2 

22.77 

94.95 

2 

93.738 

0.882 

POM 3 

3.45 

98.40 

3 

194.68 

0.915 

POM 4 

0.77 

99.17 

4 

347.35 

0.843 

POM 5 

0.55 

99.73 

6 

672.27 

0.688 

POM 6 

0.22 

99.95 

5 

512.54 

0.731 

POM 7 

0.04 

99.99 

7 

922.03 

0.785 

Companion 1 

n/a 

n/a 

12 

3623.6 

0.529 

Companion 2 

n/a 

n/a 

20 

6869.5 

0.550 

Companion 3 

n/a 

n/a 

12 

3623.6 

0.939 

Companion 4 

n/a 

n/a 

12 

3623.6 

0.792 

Companion 5 

n/a 

n/a 

12 

3623.6 

0.755 

Companion 6 

n/a 

n/a 

12 

3623.6 

0.737 

Companion 7 

n/a 

n/a 

12,20 

3623.6, 6869.5 

0.516,0.585 




Figure 34: Comparison of transverse drive point Figure 35: Close-up of transverse drive point 

displacement PSD for the 12. 8N RMS loading. displacement PSD for the 12. 8N RMS loading. 




Figure 36: Comparison of in-plane drive point Figure 37: Close-up of in-plane drive point 

displacement PSD for the 12. 8N RMS loading. displacement PSD for the 12. 8N RMS loading. 
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VI. Concluding Remarks 

An approach has been developed to utilize experimental acceleration data to guide the modal basis selection for a 
nonlinear reduced-order simulation. Two variations were explored; one in which only the transverse acceleration 
data was used, and one in which both transverse and in-plane acceleration data were used. In the former, a hybrid of 
transverse acceleration derived POMs and in-plane companions were used to select the modal basis, while in the 
latter transverse and in-plane acceleration derived POMs were used. Both variations were able to identify bases 
which resulted in favorable reduced-order simulations relative to the simulated experiment, and relative to reduced- 
order simulations whose bases were determined from the earlier approach using the simulated displacement 
response in the system identification process. Results using both transverse and in-plane accelerations in the basis 
selection yielded higher quality simulations than the hybrid approach because a greater number of in-plane modes 
were selected through the system identification process. Use of modal expansion instead of MAC for in-plane basis 
selection could benefit both in-plane POMs and companions as both were distorted relative to the displacement 
derived POMs. 

The approach was demonstrated for an idealized case in which all modeling parameters were known, and a more 
realistic case in which modeling uncertainties associated with boundary conditions and added mass were included. It 
was found that these modeling uncertainties did not negatively affect the reduced-order model results any more than 
the full-order model from which it was derived. Through this work, the groundwork has been established to permit 
the use of actual experimental data for basis selection in future studies. 
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